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Motivated by a mean field-like resistive behavior in magnetic fields commonly seen in var- 
ious superconducting cuprates and organics with strong fiuctuation, superconducting (SC) 
quantum fluctuation effects on resistive behaviors are reexamined by putting emphasis on 
their roles in the so-called thermal vortex liquid regime. By incorporating the quantum fluc- 
tuations and a vortex pinning effect in the Ginzburg-Landau (GL) fluctuation theory, it is 
found that the resistivity /5(r)-curve sharply drops, with no fan-shaped broadening, at a 
vortex-glass transition point far below an apparent upper critical field H*2{T) as a result of 
a quantum fiuctuation enhanced by an adequately small condensation energy or by a strong 
field. Fittings to p-T data of cuprates and organics are performed by phenomenologically 
including a SC pseudogap region created by high energy SC fluctuations and possible fiuctu- 
ations of competing non-SC orders. By examining La2-2;Sr2,Cu04 data over a broad doping 
range, we obtain such conclusions, consistent with recent experimental results, that the in- 
plane coherence length of hole-doped cuprates decreases with approaching the underdoped 
limit even in the presence of fluctuating competing orders and that the condensation energy 
density {Hc{0))^ = 0.5[(/>o/(27rA(0)^o)]^ is maximal near the optimal doping. Further, the case 
of disordered quasi 2D films showing the field-tuned superconductor-insulator transition is 
also examined for comparison and discussed in relation to data reported recently. 
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1. Introduction 

Macroscopic behaviors in high temperature cuprate superconductors in nonzero magnetic 
fields {H 7^ 0) have led to a renewal of our knowledge on the superconducting (SC) fluctua- 
tion. In magnetic fields typically of tesla range applied perpendicular to SC layers, both the 
resistivity p(T, H) and heat capacity in some optimally-doped cuprates show a field-induced 
fan-shaped broadening behavior near and below an apparent upper critical field H*2{T) to 
be estimated from thermodynamic quantities. Due to this consistency between p and ther- 
modynamic quantities, the resistive broadening was explained as a generic phenomenon in 
a disordered state (thermal vortex liquid regime) created by the thermal SC fluctuation in 
nonzero fields. ^"'^^ After that, the vortex lattice melting transition and its effect on the van- 
ishing of Ohmic resistivity have been studied extensively^"^-* as key issues on the vortex phase 
diagram in real systems. 
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However, researches on the vortex phase diagram of cuprates were hmited to the case with 
low enough H/H*2{0) values, in which the fluctuation is purely thermal and the vortex pinning 
efi^ect may be weaker so that the discontinuous nature of the melting transition in pure case 
may remain intact. The corresponding measurement in higher H/H*2{0) values has been real- 
ized in overdopcd cuprates under several teslas where effects of the pure vortex-solidiflcation 
are rarely seen. These materials seem to have a longer T = in-plane coherence length 
and their p-T curves have shown an apparently mean-fleld like sharp drop^'^^ particularly in 
higher fields. Through a comparison with heat capacity data,^^ however, it is clear at present 
that the sharp resistivity drop has occurred much below H*2(T) (or equivalently, T*2{H)), pos- 
sibly except at low enough temperatures, suggesting an enhanced fluctuation effect creating a 
broad vortex liquid regime in those materials. The data in refs.7 and 9 may be understood by 
merely assuming a thermal SC fluctuation and noting that the normal conductivity (J„ in the 
overdoped Tl-compounds^^ is of the order lO"^ (RqcI)^^: This value is much larger than a typ- 
ical one of the vortex flow conductivity, i.e., the mean field expression of the superconducting 
part 

as = a-an (1) 

of the total conductivity a, so that cr 2± cr„ (the normal part of a) even much below H*2{T), 
where Rq = Trh/2e^ = 6.45 (kO) is the quantum resistance, and d the distance between the 
superconducting layers in a layered superconductor. However, such a sharp drop of resistivity 
much below H*2{T) accompanied by a diminishing broadening with increasing field has been 
also observed in other cuprate superconductors^' g^^^j organic superconductors^^' with 
low H*2{T) but with much lower (T„ values (< 10{RQd)~^) and remains unexplained. Equation 
(1) implies that cr^ itself is extremely small in the thermal vortex liquid region in the tesla 
range of these materials. 

In this paper, a theory is presented to comprehensively understand such anomalous resis- 
tive behaviors in cuprate and organic superconductors with lower H*2{0). It is an extension 
of the previous work^'^) to the case with low condensation energy, in which the quantum SC 
fluctuation is not negligible in Us in 7^ but induces a resistivity curve following the normal 
(or, quasiparticle) resistivity curve Pn{T) = even below H*2{T). One of our purposes in 
this paper is to explain how to evaluate intrinsic material parameters of SC materials with 
strong fluctuation through flttings to resistivity data. Fitting results to data^^^ of underdoped 
La2-a;Srj;Cu04 (LSCO) based on the present theory were reported in ref.20. 

Here, let us briefly explain why the quantum SC fluctuation in the case is important even 
at high temperatures. To do this, it should be flrst stressed that, in a fixed H, a longer ^0 
does not necessarily imply a weaker fluctuation eflFect if recalling the fact that the Ginzburg- 
Landau (GL) fluctuation strength g2{H) in two dimensional (2D) systems near the zero field 
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transition temperature Tc is given by 



.1,6) 



« T,{\\Q)eo)H, 



(2) 



where (pQ is the flux quantum, A(0) is the magnetic penetration depth in T ^ hmit (defined 
by extraporating from the GL region), and Hq = (p^/ {2tt£^q) is the mean field upper critical 
field at T = defined at the microscopic level. A relation between i^o and H*2{0) will be 
given in the following sections. According to eq.(2), an increase of under fixed values of 
other material parameters suggests an enhanced fluctuation at a fixed H ^ 0. Note that, in 
the ordinary 2D GL theory, the Ginzburg number in H = 0^^ is independent of ^^^^ given 
by Tc[A(0)]^/(^Q(i), except a constant prcfactor, which in cuprates decreases monotonically 
with overdoping. As emphasized elsewhere, ^^-^ an increase of 52 leads to an enhancement not 
only of the thermal fluctuation but of the quantum one. In contrast to the H = case, 
the SC fluctuation in H ^ remains massive^^^ even deep in the vortex liquid region, and 
thus, its quantum contribution may play an essential role there. In Fig.l, roles of quantum 
SC fluctuation in the resistivity are illustrated (Details of calculation leading to Fig.l will be 
explained later). Figure 1 (a) also includes a comparison with optimally-doped YBCO data.^^^ 
If the quantum, SC fluctuation is taken into account in addition to the thermal one, as the 
curves in Fig.l show, an increase of §2 (in this case, of A(0)) results in a more sharp drop of p-T 
curves just above a 3D vortex-glass (VG) transition^' lying much below T*2 (-?/). Throughout 
this paper, a filled circle on each p(T)-curve indicates T*2{H). The sharp drop of p in Fig.l 
(b) is a contrast to the fan-shaped broadening around H*2 in the case Fig.l (a) dominated by 
the thermal SC fiuctuation^'^^ and is a consequence of a combination of a (pinning- induced) 
3D vortex glass transition and a reduction of brought by the quantum fiuctuations in the 
thermal vortex liquid regime. 

It has been often argued so far^^ in relation to the fluctuation effects in cuprates that the 
SC fluctuation in bulk materials with a large Ginzburg number near and in i7 7^ will be 
well described as the thermal fluctuation in the 3D XY model. Based on the above-mentioned 
fact, this conventional picture is invalid. The SC fluctuation in the limit of strong fluctuation 
is dominated in 7^ by its quantum fluctuation contribution. Then, the resistance follows 
the normal resistance curve even below H*2{T) without the fan-shaped broadening and shows 
a mean field-like sudden drop at a VG transition induced by a vortex pinning much below 



The above-mentioned reduction of ag in the thermal vortex liquid regime is due to the 



However, in the temperature range where a pinning-induced VG fiuctuation effect is negligible, 
a nearly classical (vortex flow-like) behavior intervenes these two quantum behaviors so that 



H*2iT). 



same origin as the insulating crs(T) in the quantum vortex liquid regime near T = 0.^^'^^^ 
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Fig. 1. (a) Resistivity curves (solid lines) in 4, 8, and 12(T) calculated in terms of A(0) = 0.11(/ini) 
and 12 (T) data (symbols) in twinned optimally- doped YBCO.^^^ See §3 for the details of calcu- 
lations and other parameter values used here, (b) Same curves as those in (a) except the use of 
A(0) = 0.35(/um). Each darked circle on a p(T') curve denotes T*2{H) in each H. 



they can be conveniently seen as independent ones. In fact, the latter, i.e., a reduction of 
as near T = 0, is limited in most cases to a very low temperature {T/Tc ^ 1) window 
(see Fig. 11 below) and is essential to understanding the field-tuned superconductor-insulator 
transition (FSIT) behavior^^^ which cannot occur without quantum SC fluctuation. Further, 
thermodynamic quantities such as the magnetization at such low T rapidly vary near H*2 (0) 
with sweeping H, reflecting a rise^^^ of dimensionality of fluctuation on approaching T = 0. In 
contrast, the reduction of ag, i.e., the flatterning of resistance, in the thermal regime appears, 
as the example of Fig.l shows, even at high temperatures comparable with Tc in systems 
with moderately strong fluctuation, and the corresponding T-dependence of thermodynamic 
quantities is broadened as the quantum fluctuation is stronger (see Figs. 8 and 9 below). 

Roles of quantum fluctuation of vortex positions in a perfectly clean vortex solid have 
been examined by Blatter and Ivlev^^^ as an explanation of high field behaviors of the first 
order melting line of a clean vortex solid in optimally-doped YBCO. Note that the fluctuation 
of vortex positions is included in SC fluctuations because the vortex positions are nodes (zero 
points) of the SC order parameter. However, it is well understood that the quantum effect on 
the melting transition of the optimally-doped samples is usually negligible. For instance, an 
observed field-induced deviation (reduction) of the melting line from the low field behavior is 
quite small and can be understood rather as a consequence of pinning disorder which is more 
effective with increasing fields.^^ On the other hand, the quantum fluctuation effect in the 
thermal vortex liquid regime has not been examined there. Since the fluctuation with lower 
energy becomes more important upon cooling in the thermal vortex liquid regime (i.e., the 
disordered non SC phase), it is clear that the quantum fluctuation is much more important, 
e.g., near H*2{T) rather than near the melting line. As Fig.l (a) shows, however, the quantum 
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effect on the resistivity curves, and hence, on the melting transition line, of optimally-doped 
YBCO is quite negligible. When examining in this paper resistivity curves suggestive of a 
remarkable quantum fluctuation effect, a corresponding quantum effect on a VG transition 
line replacing the quantum melting line^^^ in clean limit is relatively negligible and will not 
be examined theoretically. 

This paper is organized as follows. First, in §2, a semi-quantitative theory describing 
transport phenomena in the thermal and quantum vortex liquid regimes of 2D-like systems is 
given. Results in §2 are used in the ensuing sections by incorporating microscopic details to 
examine data on superconducting cuprates^*^^ (§3) and organics^^^ (§4), and it is shown that the 
rapid vanishing of resistivity, often seen in these materials, is a consequence of a competition 
between a strong quantum SC fluctuation and a 3D VG fluctuation. For comparison, resistivity 
curves in disordered quasi 2D systems with s-wave pairing are also discussed in §5, and it is 
emphasized through Fig. 13 that the rapid vanishing of resistivity much below H*2(T) cannot 
be peculiar to the strongly-correlated electron systems. In addition to data analysis in nonzero 
fields, effects of the SC pseudogap width and the quantum fluctuation on the critical region 
near Tc in H = are also considered in §6 and an appendix. These two ingredients are 
expected to be the origins of the unexpectedly^^-* narrow H = critical region in underdoped 
cuprates. 

In contrast to the dc electric conductivities and the Nernst coefficient, any SC fluctuation 
effect in a nonsuperconducting (non-SC) phase does not appear, at least in the GL approach, 
in the thermal conductivity.^''') The latter type of quantities are dominated by a quasiparticle 
transport^^' even much below i^o- In this paper, we focus primarily on the SC contributions of 
the former ones and assume their quasiparticle contributions (such as to be estimated from 
experimental data. Since a strong quantum fluctuation deviates as from its mean field vortex 
flow expression o"mF) we do not expect a possible change of cjmf arising from the gap nodes 
in unconventional superconductors to lead to a serious discrepancy in our results. Actually, 
the resistivity defined from the microwave surface impedance data^^^ is usually remarkably 
different from dc resistivity data and is rather comparable with ctmf derived on a single vortex 
level. 

2. Expression of Transport Quantities 

In this section, theoretical expressions useful in examining and fitting experimental data 
of the resistivity p and the Nernst coefficient N (or, the transport entropy s^) are derived. 
Although the basic framework of the theory is essentially the same as that given in ref.23, we 
need here to reexpress it in a form applicable to superconductors other than the s-wave dirty 
films considered there. Our method of renormalization of the SC fluctuation is essentially an 
extension of the Hartree approximation^' to the case with quantum fluctuation and a vortex 
pinning effect. Readers who are not interested in the details of derivation of the theoretical 



5/?? 



J. Phys. Soc. Jpn. 



Full Paper 



results may skip the main part of this section and jump to the final paragraph of this section. 
We start with a 2D GL action 



as a model of a 2D-like layered superconductor under a field H perpendicular to the layers. 
Here, Q = —iV + {2TT/(po)A(r) is the 2D gauge-invariant gradient, iP{t) = Xlcj V'<^^~*'^^ is a 
single componet pair- field (SC order parameter), (3 the inverse temperature, u = 27rn//3 with 
integer n, t the imaginary time, and b > 0. The random potential ^(r) has zero mean and 
satisfies u{r)u{r') = bp{r — r'). Although the nonlocality of bp{r) is not negligible in T ^ 
limit, we assume that bp{r) can be replaced hereafter by 6{r) multiplied by a coefficient bp 
because no fluctuation effects close to zero temperature arc considered in this paper. Further, 
the 3D nature due to a coupling between the SC layers was neglected in writing eq.(3) by 
assuming a strong anisotropy. This assumption is based on the fact found through the previous 
data fittings^' that, as far as the vortex pinning is ineffective, the effect of the layer coupling 
on the in-plane electric conductivities in 2D-like systems is extremely weak even deep in the 
liquid regime. Hence, the 3D nature will be incorporated later in considering a VG contribution 
to the conductivity. An additional dynamical term i^' uj[ipuj]* ipuj leading to a fluctuation Hall 
effect and resulting from a particle- hole assymmetry was neglected in eq.(3). This is justifled 
as far as |7'| <C 7. 

When the GL approach is applied to the low T and high H region in which any phase- 
only approach is inapplicable, -dependences of the coefficients 7, ^, and b need to be taken 
into account since the familiar low T divergences^^^ of these coefficients in zero field are cut 
off by the orbital depairing effect of the magnetic field. Thus, their Q-dependences should 
not be treated perturbatively. Although, strictly speaking, the Pauli paramagnetic depairing 
effect will also suppress the low T divergences, it will be assumed that it becomes important 
only in much higher fields than the field range which is focused on in the present work. 
Then, by expanding in terms of the Landau levels (LLs), ip{r) = ^„ '/?n.p^tn,p(i'), the 
coefficients 7 and fi of quadratic terms of eq.(3) are represented as n-dependent ones, fin 
and 7„, where n (> 0) is the LL index, p is a quantum number measuring the degeneracy in 
each LL, and n„p is an eigenfunction in n-th LL. For instance, the bare fluctuation propagator 
=< \ipn^piuj)\'^ > valid in the 6, 6p ^ limit is expressed by ^i'^^d a; |) = (7„|a;|+;[i„)^-^, 
where the triangular bracket denotes the ensemble average on The microscopic mean field 
transition point Tq{H) (or Ho{T)) is defined by fj,Q = 0. Hereafter, To(0) is often written 
merely as Tq. Detailed forms of arid /i„ will be given separately in the following sections. 
On the other hand, no n-dependence of the coefficients b and bp need to be specified throughout 
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Fig. 2. Diagrams expressing the four terms of r.h.s. of eq.(4) ((a) to (d)) and the VG susceptibility 
(f). The interaction vertex correction (the hatched rectangle) is expressed in the diagram (e), and 
the diagram (g) defines bp. The solid lines, the solid dots, the chain line, and the thin dashed line 
in (g) express, respectively, Qo, b, Qn {n > 1), and bp. 



this paper, since a high H approximation is invoked below in obtaining expressions useful for 
analyzing experimental data, and hence, the mutual interaction between the SC fluctuations 
and the interaction between a fluctuation and the random potential u are considered primarily 
within the lowest LL (LLL) with n = 0. In particular, in examining the data of cuprates and 
organics in which microscopic descriptions are still controversial, it is appropriate to assume 
that b is one of material parameters related to the magnetic penetration depth A(0), while the 
pinning strength bp is an extrinsic parameter dependent on real samples used in experiments. 

To renormalize the ^-fluctuation with low energy, the LLL approximation will be used. 
That is, a relatively high field range is assumed, in which the renormalized energy scales 
(masses) of higher (n > 1) LL fluctuations are much larger than that of LLL and do not 
deviate remarkably from their bare ones fin less sensitive to T compared with //q- Then, 
the renormalization of higher LL modes can be assumed to have already been accomplished 
independently. Further, the vortex pinning effect in the LLL-fiuctuation renormalization will 
be incorporated to the one loop order.^^^ Then, the renormalized LLL fluctuation propagator 
^o(|^<^|) is determined by 

(e^o(|w|))~^-7o|w|-Mo = So + ASi (4) 

Diagrams for the four terms of the r.h.s. of eq.(4) are given by Fig.2 (a) to (d). Here, the 
factor bp is a renormalized pinning vertex, sketched in Fig.2 (g), and is given by 

where Uy = b(3~^ Qq{\oj\) / {2'iTr'^d) . Consistently with this approximation, the VG transi- 
tioin point defined in the Gaussian approximation of VG fluctuation (i.e., the mean field VG 
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transition point) is determined according to the ladder diagram Fig. 2 (f) expressing the VG 
susceptibility^^-* as the limit Hvg,o ~* +0, where 

Although, strictly speaking, this pinning-renormalization replacing bp by hp is merely valid 
far above the vortex-solidification line in the pinning-free limit^) and is diagrammatically 
consistent just with the pinning-free fluctuation renormalization in the case with no AS; 
(see below regarding AS;), this approximation will be used hereafter for practical purposes 
because all resistivity data we will examine below belong to the cases where the first order 
vortex solid-liquid transition was destroyed by the pinning. 

The main roles of LLL mass renormalization are played by the Hartree term Sq of the self 
energy corresponding to Fig. 2 (a) and expressed as 

For the coefficient 6, the ordinary GL expression^ 

167r2A(O)V(0o^^o) in low H/Hq wih be used 
near and below T*2{H). The validity of this identification will be discussed in §6 in relation 
to the cuprates. Regarding the additional renormalization (correction) term AS; within the 
LLL, an approximation of RPA type sketched in Fig. 2 (b) 

AE,^r'E5.(H)5^iii^^§»^, (8) 

^ 2i^00(|'^|j 

(see eq.(2.11) of ref.21) will be useful below, where £'oo(|^^|) = fi~^ Qo{\uj\)Qi:){\uj + VI\). The 
term AS; is negligible in the thermal 2D case^^ and, as far as a qualitative study of resistivity 
curves is concerned, may be negligible even in the quantum 2D case. However, to make sure, 
this term will be included when attempting to fit resistivity data in §3 and 4. 

The term AS/j implies higher LL contributions to the LLL mass renormalization. It should 
be remarkcd^^^ that the higher LL fluctuations are not negligible even in the present high H 
approximation but do contribute to a downward shift of Tq{H) to T*2{H) corresponding to the 
apparent upper critical field H*2(T) , just like in zero field case^^-* where a reduction ATc of 
the mean field transition point is ascribed to the amplitude-dominated fluctuation with high 
energy. This downward shift is conveniently represented, as in Fig.2 (c), by a sum of Hartree 
diagrams. The resulting temperature Tq — ATc corresponds not to the true SC transition point 
Tc at which the linear resistance vanishes but, roughly speaking, to an onset of a sharp resistive 
vanishing or of the critical region accompanying the transition at T^. Rather than Tc, this 
resistive onset, denoted as Tco hereafter, appears in the ensuing expressions for H ^ 0. Since 
Fig.2 (c) is ultraviolet divergent even in 2D case when the quantum fluctuation is included, 
its i7-dependence is a correction, and hence, AS^ may be identified with its expression in 
H = case. Further, when the presence of other non-SC order parameter fluctuations coupling 
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Fig. 3. Schematic H-T diagram in which various crossover or transition curves defined in the text 
are drawn. The low H limit of _ffvg(r)-line corresponds to Tc- 



to the SC fluctuation is not negligible as in the underdoped cuprates (§3), they should also 
be incorporated into AS/^. It will be argued in §3 that even such non-SC fluctuations do 
not bring additional -ff-dependences. Then, AS/^ may be written as ln(To/Tco),^^'^^) and the 
r*2(-ff)-line, reducing to T^q in i/ — > limit, is determined by /^o + AS/^ = 0. In Fig.3, the 
characteristic temperatures and fields defined above are sketched in the H-T diagram. 

Due to the presence of the pinning term oc hp, solving eq.(4) selfconsistently requires 
not a difficult but a very cumbersome numerical-integration even if neglecting AS;. A more 
cumbersome situation is encountered^^-* when trying to present a practical expression on 
the VG contribution cjvg to the conductivity. Since providing theoretical formulas useful in 
analyzing experimental data is the main purpose in this section, we will not try to give a 
selfconsistent solution of eq.(4) but invoke the following approximation to give practically 
convenient expressions: The form ^o(l'-^l) = (ToI'^I + (^o(O))^^ will be assumed. That is, 
any pinning- induced renormalization of frequency dependences will be neglected. Then, eq.(4) 
is replaced by 

(^0(0))-' = ^i^ + Hn/no) + So + ASi - -^^0(0). (9) 

Consistently with this, Sq is represented as 

= ^^4^ I ^ e^+(ao(o))-- 

where the cutoff €c is a constant of order unity, and Eoq{uj) becomes 
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EooiM) = go{\io\ + \n\) go{\u;\) (11) 

1 + ^^0(1^^1/2) 

after simply arranging the a;-summation. Although, strictly speaking, this simplification is 
valid when 1 — /ig o ^ !> its justification should be discussed rather through a fitting to 
experimental data. Through computations leading to the results in the ensuing sections, we 
verified that refining this approximation was unnecessary except in s-wave disordered films at 
low enough temperatures (see §5). 

Next, bearing eq.(9) describing the SC fluctuation in LLL in mind, expressions of transport 
quantities will be derived in terms of Kubo formula for them. The superconducting part cr^ 
of the electric (diagonal) conductivity is written in the form 



d 



dcrs= { ) {je, xi-^^) je, x(i^)) 



(12) 



where the overbar (bracket ( )) implies the random (ip) average. The corresponding expression 
of the transport entropy is naively assumed to take the form^^^ 



ds^ = P(Po(^--^^ {je, x(-if^) jh, y{^^)) 



(13) 

(see, however, the next paragraph). If using the real time representation, the spatially-averaged 
electric and heat currents for the ordinary phenomenological GL model are given by 

27]- 

je,x = il — {'4^*Qxi^ + C.c)^p. (14) 
00 

3h,y = Co(i[^t^*(0]Q2/^(0 +c- C)sp, 

respectively, where t is the real time, and the bracket ( )sp denotes the spatial average. Within 
the high H approximation in which je consists of LLL and the next lowest LL, the prefactor 
^o(27r/0o) of je in eq.(14) generally becomes 7rr^/(0o ^i(O)) in the vortex liquid region irre- 
spective of the microscopic details. This fact proved in ref.23 (see eq.(2.23) there) is a direct 
consequence of gauge-invar iance. 

In contrast, the prefactor of ]h may be, at least at low T/Tq, strongly affected by a 
microscopic mechanism independent of the fluctuation property, and the correct should 
vanish in T/Tq limit. This issue in the s-wave case can be seen in ref.35. Since a correct 
magnitude of the prefactor of is unimportant for the purpose in this paper, for convenience, 
eq.(13) is used hereafter to obtain the transport energy ?7<^. However, some comments on this 
point will be necessary because, strictly speaking, a direct use itself of eq.(13) in obtaining a 
correct is not justified. Based on microscopic and thermodynamic results,^^) the factor of 
jh^y in eq.(14) has to be doubled near Tco, and a contribution —(f)oM from the magnetization 
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current has to be subtracted from the resulting Ts^, where M (< 0) is the SC magnetization. 
As is well known, however, M in high fields, where the LLL approximation is useful, is always 
given, except a prefactor, by the mean squared pair-field < \ip\'^ >= So/6, while we will 
explain below that the same thing holds true for T.s^ given by eq.(13). Consequently, as far 
as the transport energy = Ts^p in LLL is concerned, defined by eq.(13) with eq.(14) 
coincides with a microscopically valid one. 

Derivation of the term in as with no pinning-induced vertex correction, denoted hereafter 
as (Jfi, is almost the same as that in the pinning-free case.^'^^^ As already mentioned, the 
higher LLs are assumed to be inert in the LLL mass-renormalization by invoking a situation 
in high H or deep in the vortex liquid regime. In fact, reflecting the equivalence between the 
next lowest LL mode and the compressional elastic mode of the vortex liquid, {Qi{Q))^'^ has 
to reduce in the vortex liquid regime to /ii — /io insensitive to T. Further, within the high H 
approximation, fluctuation vertex corrections accompanied by an interaction between the LLL 
and the next lowest LL modes can be neglected in the dfl-expression. Then, our calculation 
of (jfl is the same as the previous one,^^-* and one finds 



^^^^^ - 2( Jo))2/3 E 



Qo{u)g.{.)[go{u) + gg,{u)\ (g^(o))-i +5(^0(0))- 

to 

r J, a ,2 1 

(ai(o))-2 



(15) 



27r 2sinh2(/3e/2) [^le^ + (go(0))-2)(7^2 + (g,(o))-2) ' 
where g = 71/70, Gi{(jj) = (7111^1 + (^i(O))^^)^"'^, and the prefactor [Qi(0))^'^ arises from the 
above-mentioned jg-vertices. As mentioned in ref.20, however, the details of ^i(O) are not 
reflected in computed results of Ufl, in most of the T and H ranges we have examined, as a 
result of the relation 71 ^o(O)/ (70^1(0)) ^ 1. For this reason, 

{gmr^ = - HO (le) 

will be assumed hereafter both above and below T*2{H). Clearly, of eq.(15) vanishes in 
r — limit. Note that the neglect based on the high H approximation of the fluctuation vertex 
corrections in cxfl does not conflict with the inclusion of AS; in the mass renormalization. 

In the realistic case with a vortex pinning, an additional contribution cTyg to ag created 
by a pinning-induced vertex correction becomes divergent on approaching a 3D VG transition 
point Tvg^'^^) from above. Near Tyg, the contribution of quantum SC fluctuation to (Tyg is 
negligible, and a civg-expression in the thermal case 

dRqayg = Cp^^^^— y^, (17) 

where t = T/Tco, will be used in comparing our theory with experimental data by choosing 
tyg = Tyg/Tco and the prefactor Cp as being sample-dependent (i.e., pinning-dependent and, in 
the case of cuprates, doping-dependent) parameters. Although the prefactor Cp should depend 
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not only on bp and b (i.e., the pinning and fluctuation strengths in = 0) but on i?,'^^^ for 
simphcity, it wifl be assumed to be i^-independent in §3 and 4. The exponent s is known to 
depend on the dimensionahty of pinning potentials dominant in the sample we focus on, and, 
strictly speaking, it is difficult to predict an appropriate value of s through each fitting far 
above T^g. Throughout the fittings to be explained below, s = 4.0 was always assumed. 

On the other hand, in situations where the system is 2D-like in spite of the presence of 
pinning effect, a true divergence of dyg may not occur, and hence the quantum fluctuation 
contribution to o^g is not negligible. As an approximate (Tyg-expression appropriate to this 
situation, the expression derived within the Gaussian approximation in ref.23 

( MoM ^ ' ^ A_ + c:^eJ-^{Q.i^) + ao(o)) (18) 



- (ao(c^))'ln 



1 + Q,o + 2|c^l7ogo(0)[^vg.o] 
l + 2|a;|7ogo(0)[evg,o]^ 



— 1/2 

will be used below, where ^vg,o = /^vg is the dimensionless VG correlation length expressed 



in unit of rg, and jUyg = /^vg,o + v7oT^o(0)- This form of ^Uyg is an expression useful for inter- 
polation near an apparent 2D quantum VG transition field B*^^^^ and the term \/ "YoTQq (0) 
plays a role of modelling the presence of the quantum VG critical regime within the Gaussian 
approximation. Further, the constant Cc of order unity is related to an upper cutoff of the 
wavevector integrals and will be hereafter chosen as Cc = 1. 

In general, a^g also vanishes inT ^ limit in i? > B*^, where B*^ is defined by /ivg,o(-B = 
B*g,r ^ 0) 0. Further, Cvg,o(5 = -B*g,r) oc T'^'^, and dRqaygiB = B*^) is approximated 
by a nonuniversal'^^'^'^^ constant at low enough T where ^vg,o S> 1. Note that eq.(18) is an 
expression valid within the Gaussian approximation and hence, may diverge m. B < B*^ at 
a finite temperature like in 3D case. As is seen later in Fig.l2, however, it is possible that it 
remains nondivergent at nonzero T, depending on the microscopic details. In the fittings, we 
will use either eqs.(16) or (18), depending upon the situations. 

One might wonder if the fact that both of eqs.(15) and (18) in 5 > B*^ vanish in low T 
limit is not a result of the neglect of some vertex corrections in the Kubo formulas. However, 
it was proved^^^ at least in the pinning-free case that all terms including the vertex corrections 
vanish in low T limit. It is trivially performed to extend the proof in ref.21 to the case with 
vortex-pinning as follows. First, as in the electron systems,^^) as far as the conductivity prior to 
the random average is considered, the fluctuation propagators t/„(|a;|) (n > 0) depend on two 
coordinates and can be represented in a form like ^n(|a;|; ri, r2) = X^^u ^^/Lt(i'i)^n,;tt(|^|)'"^(i'2) 
where u^j, is an eigen function specified by a quantum number /x. Since the proof in ref.21 
is applied in the same way as far as the spectral form (i.e., frequency dependence) in the 
low frequency limit remains dissipative and is valid irrespective of the forms of coordinate or 
wavevector dependences of Qn^ it is concluded that the sum cTfl + CTyg vanishes in low T limit 
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even if the vertex corrections are included. Therefore, by combining the high H approximation 
with this, the neglect of the vertex correction in eqs.(15) and (18) is safely valid. 

Just as for the conductivity, the transport energy [/^ may also be examined by neglecting 
the vertex corrections in the Kubo formula because the heat current is also accompanied by 
the next lowest LL mode. On the other hand, the pinning-induced vertex correction related 
to the VG fluctuation may be neglected in assuming a weak pinning because no divergent 
contribution near T^g will arise in this quantity as a result of the additional time-derivative 
in the j/^-exprcssion compared to the jg-one. By arranging the frequency summation to take 
the 0-derivative, we obtain 



dHo{o)gi{o) ^ 



+ 



1 



707i|^^l([7i^i(0)]-' - [706^0 (0)]-^)(e;o(c^) 6^1(0;))' (19) 
([71^1(0)]-^ - [7oao(0)]-')ao(a;)ai(a;) 



(7oao(0))-i + (7iei(0))-i 

- 20|a;|(go(a;))2 + ^|u;|(gi(cc;))2 
71 70 

In the present high H approximation, eq.(19) is simplified, by neglecting terms of higher 
order in 71 ^i(O)/ (70^0(0)), as 

where, as mentioned earlier, the prefactor (^^(O))"-'^ is carried by the jg-vertex. Namely, in 
the present high H approximation, in the GL region is proportional to the fluctuation 
entropy (i.e., the mean-squared pair-field) even in the quantum case, and the mean field result 
^fiiy — T/rco)/(167r^A^(0)) is expected to be recovered deep in the vortex liquid regime if 
the fluctuation has calmed down there. 

In the following sections, the above expressions of crfl, CTyg, and IJ^ are used together with 
eqs.(5), (8), (9), (10), (11), and (16) to examine experimental data of transport quantities, 
primarily, of the resistivity 

P = (C^n + C^fl + 0-vg)"S (21) 

where the normal conductivity o"„(T) should be estimated through experimental /9(T) data in 
T > Tco by, as usual, neglecting its i7-dcpendencc and is implicitly assumed to include other 
fluctuation conductivity terms excluded from the GL description. The Nernst coefficient A'^, 
measured in refs.16 and 17, is defined as 

iV=^. (22) 

<Po 

We note that, since IJs is proportional to —<\)qM in the GL region, the magnetization M may be 
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chosen in place of U(f, as a quantity to be compared with p. Even if taking account of electronic 
details of materials of interest, A(0), Tq, Hq, and T^o are independent SC material parameters 
in an ordinary microscopic description of superconductors. They appear in the coefficients of 
the GL action such as /iq and 70. On the other hand, crn(T) is highly sensitive to the sample 
purity and can be seen as being independent of the four SC parameters mentioned above. In 
order to mimic a vanishing behavior of p(T) far below Tco, a VG term CTyg, i-e., a (sample- 
dependent) vortex pinning effect on the conductivity, needs to be incorporated. If the 3D 
form, eq. (17), is used, its prefactor Cp and a form of transition line t^g{H) = T^g(H) /Tco are 
chosen to optimize the fittings. In general, the vortex pinning in real systems may occur due 
to crystal defects other than a microscopic impurity affecting an, and further, the transition 
line t^g{H) may not be described precisely in the LLL approximation assumed so far, once 
recalling a strong i7-dependence of the vortex elastic moduli and the presence in real samples 
of a small amount of other pinning sites. That is, the vortex pinning should be regarded 
as being independent of cj„ and is not necessarily described in terms of a single extrinsic 
parameter. If the 2D form (18) is a more appropriate cJvg, the pinning strength bp is the only 
extrinsic (sample-specific) parameter at nonzero temperatures. 

3. Cuprate Superconductors 

In this section, the theoretical expressions in §2 will be applied to experimental data of 
superconducting cuprates. Since one purpose of examining resistivity data of cuprates is to 
correctly understand the doping dependences of fluctuation effects and of material parameters 
of cuprates, we will primarily examine resistivity data of LSCO of which an extensive doping 
dependence can be seen in the literature. As shown in ref.20 where the expressions obtained 
in §2 were applied, data of other quantity measured consistently are also needed, together with 
resistivity data, to correctly estimate material parameters in cuprates with small condensation 
energy. In underdoped cuprates, such a set of data measured consistently are not known 
except the LSCO data in refs.l6, 17, and 39. Below, we will proceed further the analysis to 
the overdope side and comment on other cuprate materials on the basis of available data.^^'^^^ 

Let us start with incorporating microscopic ingredients into the GL description. Since 
experimental data in several teslas are primarily examined for materials with a much lower 
H*2(0) than the optimally-doped YBCO, we need a reasonable functional form of H*2(T) which 
may not be approximated in several teslas as the ordinary linear straight line. Further, the 
time scales 7„ need to be calculated consistently with this H*2 (T) . To this end, we invoke the 
ordinary clean limit^°^ in order to describe pn and 7^ consistently (see §2 on their definition). 
For simplicity, let us assume, as in the weak-coupling s-wave pairing case, a circular Fermi 
surface. Then, they are given by 




(23) 
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Ain = ln — +/ ds ? ' . , 

\ToJ Jq sinh(s) 



respectively, where Uc = Toy/H/{2H^e^/T, L„(a;) is the n-th order Laguerre polynomial, and 
7 = 0.5771 is the Euler constant which should not be confused with the time scales jn-^^''^'^^ 
It is valuable to comment on the fact that 7„ with any odd integer n approaches zero in 
T ^ 0. Since an equilibrium vortex solid state with no Pauli limiting effect is represented by 
the LLs with even n, a dissipative vortex flow motion is created by other odd LLs^^^ so that 
the vortex flow conductivity is proportional to with an odd n. Hence, this result suggestive 
of dissipation-free vortex flow at T = in clean limit may be rather expected and significantly 
enhances the quantum resistive behaviors in cuprates and organics at finite temperatures. 

As noted in §2, the apparent upper critical field line H*2{T) approaching T^o in low H 
limit (see Fig.3) is determined by //q + AS/^ = 0. When AS/^ = ln(ro/rco), we obtain 



H:,{T) = Ho(^^ym, (24) 

from eq.(23), while Ho{T) is given by Ho^{T/To), where the function satisfies $(0) = 1 
and $(1) = 0. In particular, a large enough value of the parameter Tq/Tcq significantly affects 
fluctuation phenomena in nonzero fields. 

Note that, when Tco rather than Tq is chosen as a temperature parameter scaling T, Hq 
is replaced by H*2{0) = Hq{Tcq/Tq)'^ . By combining this with eqs.(4) and (11), we find the 
property 

afl(ro/reo,//o,A(0)) = afl(l,/f*2(0),rcoA(0)/ro) (25) 

to be valid in the high H approximal ion. Namely, if the SC pseudogap region with the width 
^0 — is neglected, = V^o/(27ri?o) is overestimated, while A(0) is underestimated. A 
discussion based on this fact will be given in §6. 

As the first example of comparisons with data, the case of optimally-doped YBCO shown 
in Fig.l (a) will be explained. In this case, the p{T) curves in the tesla range show the fan- 
shaped broadening in the thermally- induced vortex liquid region below Tco, and the thermal 
SC fluctuation reduces /9-values near T*2{H) compared with the mean field result of resistivity. 
Figure 1(a) includes 12 (T) data in ref.22 and p{T) curves computed in terms of the material 
parameters given in Table I, while the curves in Fig.l (b) result from the same set of material 
parameters, except A(0) = 0.35 (/um), as in Fig.l (a). Further, the normal resistivity cr"^ = 
0.135T/rco(m0.cm) and the 3D-like o-yg, eq.(17), with Cp = 1.3 x 10"^ and tyg = 1 - h - 
1.2h^/^, where h = H/H*2{0) <^ 1, were assumed. Since, in several teslas, the /i- values of 
optimally-doped YBCO are much lower than those in most of underdoped and overdoped 
materials including LSCOs to be discussed below, an additional contribution^) composed only 
of thermal higher LL modes to af was also included in obtaining the solid curves together 
with the quantum contribution accommodated in eq.(15). In Fig.l (a), however, a change of 
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Fig. 4. ?7^-curves (solid curves) at 3, 6, 8, and 12 (T) obtained consistently with those in Fig.l (a) 
and 12 (T) data (open squares). In this and other figures of U^{T), is represented in unit 
of 10-i2(J/m). 



resistivity value brought by the addition of higher LL contributions is within several percents 
in magnitude, and the quantum contribution was quite negligible. The positions of H*2{T) for 
each resistivity curve are denoted by filled circles both in Fig.l (b) and the figures appearing 
hereafter. Further, we note that, although a small but nonvanishing Tq/T^o — 1 was taken 
into account in Fig.l favoring an agreement with the consistent data of (see Fig. 4), the 
obtained values of A(0) and remained unchanged compared with those in previous fittings^^ 
where the quantum contribution was neglected from the outset. 

We note that Tq wc define is the bare mean field SC transition point defined prior to 
including a coupling with possible non SC orderings competing with the SC ordering. We 
expect two types of origins of the SC pscudogap region Tq — T^o in cuprates. First, inclusion of 
a fluctuation of a competing non-SC ordering, such as a spin or charge ordering, coupling to 
the SC order parameter can lead to a reduction of the mean field (MF) SC transition point. 
The resulting MF transition temperature in if = will be called hereafter. This can be 
modelled by adding the term 

S^=u^ JJ d\\<Pns\^\i^\^ (26) 

in GL action, where (pns denotes a non-SC order parameter fluctuation, and only a fluctuation 
competitive (or repulsive) to ijj is assumed here through the condition Ux > 0. Since |^nsP 
is replaced by the averaged value < |0nsP > in constructing an effective action on ijj, 
will be expressed as Toexp{—Ux < |0nsP >) (< ^o); and consequently, we has only to replace 
ln(r/ro) in fj,o by ln{T/T^^), where < \(f>ns\'^ > was assumed to be if-independent. The 
second origin of To-shift is nothing but AS/j in the LLL mass renormalization outlined in §2, 
where it was assumed to arise entirely from the SC fluctuation in higher LLs. An additional 
contribution to AS^ also arises from Sx and similar higher order coupling terms between ip- 
fluctuations and (f)ns- As in §2, by assuming this contribution to AS/j of such a 0n5-fluctuation 
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Fig.l(a), 4 Fig.5 Fig.8 



A(0) (^m) 


0.11 


0.43 


1.9 


Ho (T) 


272 


235 


330 


To (K) 


120 


96 


90 


Tco (K) 


92 


32 


15 


d (nm) 


1.5 


1.5 


1.5 



Table I. Material parameter values used in Figs.l, 4, 5, and 8. Parameters related to the pinning 
effect and cr„ are given in the text. 



0.45 




Fig. 5. Fittings to LSCO x = 0.08 data at 12 and 26 (T) of (a) p and (b) N in Ref.l6. 



to be also iJ-independent, AS^ can be identified with ln(r^^^/rco). That is, T^^ does not 
appear in hq + AS^. 

The assumption of i?-independence of the ^n^-fiuctuation can be justified by comparing 
with experimental observations as follows. First, according to NMR measurements^^) in the 
pseudogap regime, the antiferromagnetic (AF) fluctuation, which is competitive with the SC 
fluctuation and should be dominant in the NMR signal in the nodal directions, is suggested 
to be insensitive to H at least in T > Tcq. On the other hand, a local AF ordering was 
shown to be enhanced with increasing H primarily below the irreversibility line where the 
ohmic resistivity is absent.^^^ This iJ-dependence of AF ordering is essentially linear in H, 
reflecting the number of vortices, and possibly, will be a consequence of the spatial variation 
of IV' I near the vortex cores.^^'^^^ However, in the high H approach of GL theory, the spatial 
variation of |V'| is reflected only in /3a — 1 which is negligible^) in discussing thermodynamics 
and transport phenomena above the irreversibility line (i.e., in the vortex liquid regime). 
Therefore, we believe that the assumption of i7-independence of ^n^-fluctuation is valid in 
the vortex liquid regime. Even if the <;!!)„s-fluctuation induces an i^-dependence in the effective 
action on the ■^-fluctuation, both fluctuations are competitive with each other, and hence, a SC 
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Fig. 6. Resistivity data (symbols)^) of LSCO in the same temperature range at three doping levels 
(a) X = 0.08, (b) 0.15; and (c) 0.2 and corresponding theoretical results (solid curves). The darked 
circle on each p{T) curve denotes T*2{H) in each H. 

coherence length defined from a iJ-dependence in the effective action would be enhanced by 
the presence of <;ii„s-fiuctuations. As is repeatedly seen hereafter, however, ^o-values resulting 
from the fittings to data of underdoped cuprates already become short enough. For this reason, 
the extra iJ-dependence due to the (?!)„s-fiuctuations is likely to be quantitatively negligible, 
and hence, we assume the contributions of ^„5-fiuctuation to have been fully accommodated 
as ln(ro/T^^). Although, as a result of this, does not appear in the r.h.s. of eqs.(4) 

and (9), the mean field transition field in the presence of <;ii„s-fluctuations will be given, 
consistently with eq.(24), by Hq{T^^ /To)'^^{T /T^^) which may decrease with underdoping. 
Note, nevertheless, that the genuine is defined through the microscopic field scale Hq 
irrespective of (?!)„5-fiuctuations. 

Based on the inclusion, explained above, of competing ^^^-fluctuations, the fitting results 
to LSCOx data with x = 0.06 and 0.08 reported in ref.20 will be commented on here. Hereafter, 
in considering the doping dependence of cuprates, the spacing d between the SC layers will 
be fixed. In ref.20, the formulas given in §2 were applied by neglecting the difference between 
To and T^^. Regarding cTyg, eq.(17) with Cp = 10~^ was used for x = 0.08 case, and eq.(18) 
with bp = (27r^Qd)0.023 was applied to x = 0.06 case . According to the results in ref.20. 
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Fig. 7. C/^-curves corresponding to p-curves in Fig. 6 (c). 



ro(x) became insensitive to x in x < 0.1 in contrast to T^,, defined in ref.26 as the onset of SC 
Nernst signal, which decreases with decreasing x (< 0.1). Here, in relation to Fig.6 (a) below, 
the fitting results to the x = 0.08 data given as Fig.2 in ref.20 will be shown again in Fig.5, 
where {(Tn)~^ = 0.02905 ln(l. 6 x lO^Tco/T) (mfi.cm) and the parameters listed in Table I were 
used. Regarding the x = 0.06 case, it has been noticed recently that, even if assuming Tq to 
increase with decreasing x just like Ti, in Bi-compounds,^^^ the data can be quantitatively 
explained. For instance, the parameter values A(0) = 2.3(/im), Tq = lOO(K), Tco = 13(K), 
Ho = 493(T), and {an)~^ = 0.09941n(186.3rco/r) (mJ^.cm) were used for the x = 0.06 case, 
and the obtained curves of p and N were almost the same as those in Fig.l of ref.20 where 
a smaller Tq value, 96(K), was assumed. Contrary to this, no choice of Tq decreasing with 
underdoping has resulted in a consistency with the Nernst data^^^ particularly in the x = 0.06 
case. What corresponds to Ti, in low H limit is the mean field transition point T^^. More 
accurately, because the Nernst signal due to Gaussian SC fluctuation is usually nonzero above 
the microscopic transition point, T^^ may lie slightly below Ti,. Hence, T,^ in x < 0.1 of LSCO, 
showing the unexpected doping dependence, will correspond to T^^ affected by competing 
non-SC fluctuations becoming stronger with underdoping. That is, we expect Tq in LSCO to, 
contrary to T^^, increase with underdoping like in the Bi-compounds.^^^ 

Next, in relation to Fig.5, the doping dependence of resistivity curves will be examined in 
order to corroborate that doping dependences of material parameters predicted in ref.20 are 
valid over a wider doping range including the overdoped side. Previously, a systematic study 
of resistive behaviors at various doping levels of LSCO was reported,^) and the data in ref.8 
will be used here. We focus on the temperature range 10 < T'(K) < 35 and approximate each 
[cr„(T)]~^ curve there via a T-linear curve. Among various data in ref.8, we have examined 
three doping levels, x = 0.08 (underdoped case), x = 0.15 (nearly optimal case), and x = 0.2 
(slightly overdoped case). Fluctuation effects in strongly overdoped materials (x > 0.24) 
are expected to be quite weak, and those data will not be considered here. First of all, the 
X = 0.08 data have been fltted by assuming a Tq/Tcq value similar to that in Fig.5. Next, the 
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X 


0.08 


0.15 


0.2 


A(0) (^m) 


0.57 


0.46 


0.417 


^^0 (T) 


245 


190 


55.1 


To (K) 


100 


86.5 


40.5 


TcQ (K) 


30 


31 


30 


a (nmj 


1.5 


1.5 


1.5 


(mO cm) 


0.245 + 0.27i 


0.19 + 0.0426t 


0.055 + 0.0312i 


Cp 


0.04 


0.01 


0.033 


ayg 


0.5 


0.67 


0.67 


/3vg 


3.4 


0.8 


2.26 


HM (T) 


0.352 


0.385 


0.23 



Table II. Parameter values used in computation of resistivity curves in each figure of Fig. 6. The 
parameter t denotes T/Tco- 

parameter values for x = 0.2 case were chosen by favoring a semiquantitative agreement with 
Nernst data in ref.l7 (see Fig.l there). Then, in x = 0.15 case, a Tq/Tcq value intermediate 
between the 0.08 and 0.2 cases was assumed. However, according to eq.(25), the key parameter 
in comparing with p data is the product A(O)^O) that is, i?c(0) = \/ 4>qHq / Anr (see 

eq.(2)) which is independent of Tq in GL theory. The obtained curves are shown together 
with the data^) in Fig.6, and the values of material parameters used for fittings are listed 
in Table.IL The parameters ayg and /3vg are included in the assumed form of VG transition 
curve tvg(^) = Tyg{h)/Tco = (1 - h)/{l + ^vg^""«), where h = H/H*2{0). This expression of 
tvg{h) is based upon the VG transition line in the mean field approximation derived within 
the LLL.^) The factor 1-h arises from the mean-squared amplitude of the pair-field, while the 
exponent ayg depends upon the dimensionality of SC fluctuation and takes 0.5 and 0.67 in its 
2D and 3D limits, respectively. This ayg is a direct measure of the sample anisotropy in the 
present approach applicable primarily to 2D-like materials (see the second paragraph of §2). 
The Qjyg values shown in Table II are compatible with the well-known fact that the hole-doped 
cuprate materials are more 2D-like with underdoping. By contrast, /3yg is sample-specific and 
also depends on both of the fluctuation strength and the pinning strength. Furthermore, there 
are at least two ingredients affecting the form of tyg{H) in realistic cases. One is an ingredient 
independent of pinning effects and induces a deviation of tyg-form from its LLL expression. 
Physically, a detail of the vortex elasticity, such as the i7-dependence of the shear modulus of 
a vortex lattice defined locally, will affect the form of tyg{H), and hence, it is not surprising 
that such a difference of the local shear modulus in ^ < 1 from its LLL expression affects the 
tyg(iy)-expression. Secondly, a small amount of line-like (or plane-like) pinning disorders in real 
systems also affect the transition line. For instance, in real systems including both point-like 
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and line-like pinning disorders, the functional form of t^g{h) depends even on the ft.- values. 
Thus, in fitting to resistivity data, it may be rather necessary to include a possible deviation of 
the tvg(i/)-form from the expression in LLL and in the purely point disorder case. To reduce 
the fitting parameters as far as possible, we have assumed that such involved ingredients will be 
incorporated by taking (3^g as a sample-specific fitting parameter independent of other material 
parameters. Nevertheless, it should be stressed that the details of the extrinsic parameters Cp 
and tvg(^) should not be important for our purpose in this section of clarifying the doping 
dependence of intrinsic material parameters which are determined, roughly speaking, through 
the (pinning-independent) upper-half of the resistivity curves. 

In lower fields in x = 0.08, the resistivity curves have the fan-shaped broadening suggestive 
of a dominance of thermal fluctuation over the quantum contribution. One might wonder if 
the result in Fig. 6 (a) that the thermal (fan-shaped) behavior in lower fields is more evident 
with underdoping is consistent with the result^") on strongly underdoped {x < 0.1) cascs^^'-^^^ 
where the thermal behavior was lost with underdoping in x < 0.1. This apparently conflicting 
result is resolved as follows. In the case of underdoped cuprates in tesla range, a strong SC 
fluctuation near Tq is weakened to some extent upon cooling down to Tco and behaves, much 
below To, like that in effectively higher fields. Roughly speaking, such a reduction of fluctuation 
arises from an increase, upon cooling in Tco < T < Tq, of the coefficient of the gradient term in 
the H = GL-expression. That is, the enhanced quantum fluctuation due to A(0) increasing 
with underdoping is partly compensated by the existence of a large SC pseudogap region 
which, in turn, results in weakening the fluctuation near and below Tco- Hence, due to a 
T-depcndcncc of an electronic origin, a long A(0) may not necessarily result in a quantum 
fluctuation-dominated behavior near Tco in the tesla range. Actually, in contrast to LSCO in 
X > 0.1, the quantum fluctuation effect on p{T) seems to be monotonously enhanced with 
underdoping in any H in the case of YBCO.^^''^^^ 

In Fig. 6 (a), no H*2{T) position (dark circle) on the 29(T) curve was indicated. Actually, 
if*2(0) for Fig.6 (a) is close to 21(T), and hence, this 29 (T) curve is an example of the case 
in which the pinning- induced drop of resistivity occurs above H*2(T) as a consequence of a 
broad SC pseudogap region. It is not surprising because a VG transition can occur anywhere 
in T < To(H), i.e., as far as the SC fluctuation is present. A similar feature will be discussed 
again in §6. 

In contrast, the resistivity curves in x = 0.2 and 0.15 cases always show a sharp drop. 
As suggested in ref.17, the high field p curves in x = 0.2 case show a sharp drop far below 
H*2(T) (dark circle) suggested from Nernst data, and hence the situation is likely to be 
similar to Fig. 1 (b) in Introduction. Actually, the much weaker iJ-dependence of p around Tco 
in X = 0.2 compared with that in x = 0.15 case is apparently inconsistent with an expected 
growth of ^0 accomanying the overdoping and rather reflects a fluctuation enhanced with 
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overdoping. This is due to A(0)^o increasing with overdoping starting from the optimal-doping, 
although the H = fluctuation is weakened with overdoping (see §1). The ^J^-curves computed 
consistently with Fig. 6 (c) are shown in Fig. 7. The obtained [/^ values and //-dependences 
semiquantitatively agree with the data shown in Fig.l of ref.l7. Properties of the [/^ curves 
are typically thermal, and, at least below 12 (T), there are no remarkable quantum fluctuation 
effect on C/^. Regarding the x = 0.15 case, we note that the VG transition position t^g{h) is 
closer to H*2{T)-\\ne compared with those in a; = 0.08 and 0.2 cases, and that, due to this, 
the sharp p-drop in a; = 0.15 case is rather similar to that in the mean field-like case with 
negligible vortex liquid region. This narrow vortex liquid regime suggested by the x = 0.15 
data may be realized if the effective pinning strength relative to the fluctuation strength is 
maximal near the optimal doping. 

One of important consequences arising from the fittings is that [i/c(0)]^^ oc A(0).^05 i-e, 
the fluctuation strength in fixed //-values (see eq.(2)), is mimimal near the optimal doping 
and, just like (Tco)^^, increases with both underdoping and overdoping from the optimal case. 
This is qualitatively consistent with the doping dependence of condensation energy density 
[//c(0)]^/(47r) estimated from the heat capacity data.^^-* It is easily understood by recalling the 
discussion on cq.(2) in §1 that the doping dependence of p{T) curves mentioned in §1, including 
the fact that the fan-shaped resistive broadening is typically seen only near the optimal doping, 
is a reflection of this doping dependence of //c(0). Further, as Table II shows, the in-plane 
coherence length to be defined from Hq monotonously decreases with underdoping over all 
doping ranges including x < 0.1.^'^-' This conclusion cannot be reached once the presence of 
the SC pscudogap region widening with underdoping is neglected.'^'') 

Examples of p and C/^ curves in a case with very large strengths of both the SC fluctuation 
and pinning effect are shown in Fig. 8. They have been computed by bearing very underdoped 
Bi-2201 data^^^ in mind (see Fig. 2 in rcf.39). A quantitative comparison with the data will not 
be attempted here because the data were taken on a film sample with a broadening of p{H = 0) 
curve over 10 (K). This sample-specific broadening at // = should be also reflected in low H 
curves of resistivity and make comparison of computed curves with the low H data difficult. 
Nevertheless, we expect semiquantitative features of the Bi-2201 data except p-curves below 
6(T) to be comparable with Fig.8. As well as in the LSCO case with x = 0.06,^°^ the 2D a^g 
expression, eq.(18), was used with bp = {27rC^d)0.15. In Fig.8, (c7„)-^ = 0.211n(280 Tco/T)) 
(mfi.cm) was used together with the parameters shown in Table I. Qualitatively, the features 
are also similar to the LSCO data in a; < 0.06.-'^^'^°^ First, the quantum SC fluctuation in Fig.8 
is strong. In fact, T*2 at 8(T) is close to 5(K) where the p(r)-curve is insulating. More notably, 
the resistivity curves above 5(K) suggest a 2D FSIT with B* ~ 6(T), like in Fig.l of ref.20. 
However, the FSIT behavior in Fig.8 is more remarkable compared to the LSCO case with 
X = 0.06.20) This is a consequence of the larger value of pinning strength bp. Namely, a large 
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enough value of pinning strength is needed together with a strong quantum SC fluctuation to 
obtain a more remarkable FSIT behavior visible even near Tcq. On the other hand, the larger 
pinning effect enhances the VG transition field at low T. Actually, H*2{0) in Fig.8 is less than 
10(T) at which a rapid drop of p at a finite T still occurs. 

Regarding the transport energy, the [/^-values in Fig.8 (b) are quite low and should be 
compared with other figures of [/^ shown in this paper. As well as the LSCO case,^°^ the 

or AT- values are two order of magnitude lower than a value expected by ignoring the SC 
pseudogap region (i.e., when Tq = Tco) and are consistent with the data.^^^ 

This comparison with underdoped Bi-2201 data corroborates the argument in ref.20 based 
on the LSCO data as follows. Since, in contrast to the LSCO case in a; < 0.1, in Bi2201 
seems to monotonously increase^^) with underdoping, one might expect underdoped Bi-data 
to behave in a qualitatively different way from the LSCO data. As seen above, however, 
the present theory in which corresponding to Ti, plays no essential roles explains the 
similarity^^' in behaviors near and below Tco of both p and C/^ data between Bi-2201 and 
LSCO. It implies that the examples of LSCO studied in ref.20 can be seen as generic behaviors 
of strongly underdoped cuprate materials below Tco- 

In this paper, the case of underdoped YBCO is not examined in details because no com- 
parable data of p and N at the same doping level in underdoped YBCO have been reported. 
Actually, the 3D nature of SC fluctuation should be incorporated in theoretical descriptions in 
contrast to other cuprates which are commonly much more 2D-like, and hence, the approach 
in §2 may not be directly applicable to a quantitative study of doping dependences of YBCO. 
Nevertheless, the following features are suggested from available p and N data: First, the fact 
that, in contrast to the LSCO x = 0.08 data^^ showing the fan-shaped broadening below 8 
(T), the p data in underdoped YBCO^^'^^'^^^ entirely show a sharp drop implies that the en- 
hancement of quantum fluctuation always overcomes a reduction of fluctuation upon cooling 
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arising from a broad SC pseudogap regime, and hence, a Tq — Tco value much smaUer than that 
in LSCO is expected. This narrower SC pseudogap regime is presumably consistent with the 
Nernst coefficient^^^ larger than in LSCO and Bi-2201. In any case, the doping dependences 
of Tjy and Tq in YBCO, which we expect will be remarkably different from each other, should 
be clarified elsewhere. 

4. Organic Superconductors 

Previously, the resistive behaviors in the vortex liquid regime of k-(ET)2 organic super- 
conductors have been studied in parallel with those of cuprate materials. Typical data are seen 
in refs.18 and 19. Surprisingly, the resistivity curves in k-(ET)2Cu(NCS)2 with a wider vortex 
liquid regime have shown a sharp drop near the irreversibility line in all fields shown there. In 
K-(ET)2Cu[N(CN)2]Br with a narrower liquid regime, the resistivity p{T) curves have shown 
a clear i?-dependent crossover from the familiar fan-shaped broadening in lower H into a 
sharp drop in higher H near the irreversibility line lying much below H*2{T) estimated from 
the magnetization data. These features are much the same as those seen in cuprates^'^^'^^^ 
and are consequences of quantum SC fluctuation becoming more important as the fluctuation 
is stronger. Here, in addition to reproducing and discussing the fitting result, the low T 
behaviors will be examined. Numerical results in this section are also based on the use of 
eqs.(23) and (24) within the framework of §2. 

Resistivity data of K-(ET)2Cu[N(CN)2]Bri8) are fitted in terms of the present theory, and 
the results are shown in Fig. 9 (a). The 3D form of civg, eq.(17), is assumed together with t^g = 
{1 — H / Hq)"^ / {1 + 3.52{H / Hq)^^'^) and Cp = 0.0035, and the material parameters used there are 
A(0) = 0.72(/im), To ~ Tco = 12 (K), Hq = 18.3(T), d = 1.5(nm), and the normal resistance 
i?„(0)= 0.0574 + 0.0243 1^/^. As a result of the large A(0)-value, the quantum SC fluctuation 
becomes essential with increasing H and leads to high field curves following Rn{T) even in 
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Fig. 10. Data of resistivity of k-(ET)2Cu(NCS)2 near T = and in strong fields-^^) 




Fig. 11. Resistivity curves computed for a qualitative comparison with the data in Fig. 10. 



T <^ Tco- Further, the Cp-value and the tyg-form suggest a much weaker pinning effect than 
in LSCO. For comparison, the [/^ curves obtained in terms of the same set of parameters are 
shown in Fig.9 (b). Note that, within the LLL, is equivalent to the magnetization, and hence 
that Fig.9 (b) can be also regarded as a typical example of thermodynamic quantity. A crossing 
behavior just below Tco, which is familiar through magnetization data in many optimally- 
doped cuprate materials, is seen in lower fields below 2(T). It implies that the fluctuation 
property below 2(T) is purely thermal. In contrast, the C/<^(r)-curve in 6(T) deviates from 
the crossing behavior in lower fields and is anomalously broadened, refiecting the stronger 
(quantum) fiuctuation. Such an additional broadening of thermodynamic quantities due to 
the quantum fiuctuation is an opposite trend to the corresponding behavior near T = 0, where 
a broadening of such quantities diminishes refiecting a rise of dimensionality of SC fluctuation 
in the quantum regime near T = 0.^^^ This is why, as mentioned in §1, the quantum behavior 
at high temperatures should be distinguished from that near T = 0. The above features of U^j, 
seem to be consistent with the magnetization data in ref.18. 

Next, the resistive behavior near T = will be discussed based on the data of k- 
(ET)2Cu(NCS)2 shown in Fig. 10 which corresponds to Fig.3 of ref.49. In the highest field 
8.59(T) above Hq defined through a study of Shubnikov-de Haas effect,^^) p(T) shows a T- 
dependence indicative of a metallic normal state and approaches a residual value Pn{0)- In 
lower fields just below Hq, p{T) first decreases upon cooling as a result of thermal SC fluc- 
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Ru/Rq 0.25 0.4 0.45 0.57 0.625 0.67 
B*JHq 0.8 0.76 0.72 0.68 0.64 0.61 
Rc/Rn 0.98 0.98 0.92 0.83 0.75 0.65 

Table III. Ru/Rq dependences of the normalized critical resistance Rc/Rn and of the normalized 
critical field B*/Hq, estimated for each case of the calculated resistivity curves illustrated in 
Fig.l2. 

tuation, while it begins to increase on further cooling and approaches Pn(0) (The ultimate 
drop of p due to the VG fluctuation at a low T of resistance below 7(T) will not be consid- 
ered here). This insulating behavior occurs even much below Ho and hence, is a phenomenon 
of a SC origin. This feature in Fig. 10 is a direct evidence of an insulating behavior^^'^'^^ 
arising from a purely dissipative quantum SC fluctuation in cases with metallic normal re- 
sistance. For comparison, we give in Fig. 11 examples of computed p{T) curves with such a 
fluctuation-induced insulating behavior at low enough T. In Fig. 11, the 2D a^g, eq.(18), was 
used, and the parameter values, A(0) = l.l(/im), bp = (27r^Q(i) 0.06, Tq = Tco = 25(K), and 
Ho = 19(T), were chosen. Further, the relations (i?Q(io-„)"i = 0.35(1 + 0.05{T/TcoY'^'^) and 
Ho = T/Tco — 1 + H/Ho were used, and, for simplicity, the vertex correction to the pinning 
strength was neglected by setting bp = bp. 

5. s-wave Dirty Films 

In ref.23, we have proposed a theory of field-tuned superconductor-insulator (FSIT) be- 
havior in homogeneously disordered thin SC films with s-wave pairing on the basis of a familiar 
electronic model in dirty limit including effects of a repulsive mutual interaction between elec- 
trons. It has been argued there how, refiecting differences in T-dependences between various 
components in a, the resistance value Rc on an apparent FSIT field B* and the resistive 
behavior around B* are affected by the normalized value Rn/Rq of the (high temperature) 
sheet resistance. However, no detailed computation results based on the derived GL action 
were given there. Motivated by a recent finding on R^ v.s. i?„ relation,^^) some computed re- 
sults on resistance curves consistent with the experimental observations^-*^' will be presented 
here. In contrast to the previous two sections, we take account of the fact that the resistance 
data of s-wave amorphous films are conventionally discussed in terms only of the Rn/Rq value 
and will not try here to fit to real data. Actually, the expressions (see Appendix A) we use 
here for microscopic parameters were derived from the simplest extension of the BCS model 
to the case with both the disorder and a repulsive interaction between the electrons and may 
not explain quantitatively materials, for example, with a strong spin-orbit scattering. 

The coefficients of each term of the GL action for the s-wave SC films were studied 
elsewhere.^^'^^^ In Appendix A, they will be given in a form useful for numerical computations. 
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Fig. 12. Calculated resistivity curves imagining s-wave dirty films for three Rn/Rq values : (a) 
Rn/Rq = 0.4 in H/Hq = 0.87, 0.82, 0.79, 0.77, 0.76, 0.74, 0.7, 0.65, and 0.6, (b) RrjRq = 0.57 
in iJ/ffo = 0.9, 0.85, 0.8, 0.75, 0.7, 0.67, 0.64, and 0.6, and (c) Rn/Rq = 0.67 in H/Ho = 0.9, 0.8, 
0.67, 0.63, 0.61, 0.59, 0.55, and 0.51. 



Applying those GL coefficients to eq.(9) and using eq.(18) as Uvg, we liave examined tfie 
resistance curves near in details, and tfie results are shown in Fig. 12. Since, as mentioned 
in Appendix A, the value of the parameter tTq was fixed, the normalized sheet resistance 
(defined at high T) Rn/Rq is the only material parameter for our calculation of Rgda and 
determines both strengths of quantum fluctuation and vortex pinning. Let us define B* as the 
field at which R{T) takes a value Rc insensitive to T at low enough T within the examined 
temperature range. 

Prom the resistance curves shown in Fig. 12, one will find the following two features on the 
curves near B*. First, for lower values of Rn/Rq, the curves near B* show an insulating T- 
dependence at intermediate temperatures, while they, for large Rn/Rq values, rather decrease 
upon cooling in the same temperature region. This feature has been seen in various data^^~^^^ 
and justifies the scenario predicted in ref.23. At least within the present model of s-wave thin 
films, the decrease of R{B*) upon cooling for large enough i?„/i?Q-values is dominated by 
an enhancement of pinning strength, while its increase in smaller i?„-values occurs primarily 
reflecting an enhancement of quantum SC fluctuation. For instance, the insulating R{B 2± B*) 
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for Rn = 0.4i?Q is controlled primarily by the quantum behavior of (jfj rather than a^g. Next, 
as listed in Table III, Rc at lower i?„ (< 0.5i?Q) values roughly coincides with while Rc 
in Rn ^ 0.6 rather decreases with increasing The relation Rc — Rn at low disorder (lower 
Rn) coincides with the experimental results in i?„ < 5(kr2) summarized in Fig. 4 of ref.52, 
while the saturation and reduction of Rc at higher disorder qualitatively agree with Fig.4 of 
ref.51. 

Finally, we will comment on the validity of eq.(9) used in obtaining Fig. 12 on the basis of 
the microscopic informations in Appendix A. When replacing eq.(4) with eq.(9), the pinning 
effect was underestimated. This simplification needs to be reconsidered in a close vicinity 
of a VG transition point (i.e., ^vg +0). In trying to fit various data in §3 and 4, we 
did not have to discuss the details of resistivity data in the vicinity of the VG transition. 
In contrast, describing in details the resistive behaviors in 2D disordered thin films at low 
enough T and near B* (i.e., near the quantum VG transition) is needed to clarify the physics 
of FSIT. We have partly carried out numerical computations based on eq.(4), although such a 
cumbersome analysis will not be presented here, and found that, roughly speaking, numerical 
results following from eq.(9) in the cases of i?„/i2Q-values used in Fig. 12 are quantitatively 
unreliable only in T/Tco < 0.05. For this reason, no data in T/Tco < 0.05 have been shown in 
Fig. 12. 

6. Comments and Discussion 

First, a view extended over a wider temperature range of the resistivity curves in Fig.l2 
(b) is shown in Fig. 13. From this figure, a T = 2D VG transition point B* is suggested to 
apparently lie much below H*2{0), because B* itself is lowered by the quantum SC fluctua- 
tion,^^) while the FSIT behavior appearing around B* is a reflection of the VG fluctuation. 
Since the flat or insulating resistive behaviors appear in H > B*, the temperature at which 
R{T) in a field slightly below B* drops inevitably lies far below T*2{H) indicated by a dark cir- 
cle in the figure. Since the FSIT behavior is a consequence of strong quantum SC fluctuation, 
this example means that the sharp resistive drop much below H*2{T) in systems with strong 
quantum fluctuation is not an artifact of approximations used in calculations but a generic 
feature occurring commonly in clean and dirty limits irrespective of calculation methods of 
R{T). 

As shown in §3, the quantum SC fluctuation plays important roles in many cuprate ma- 
terials under a high field. If the quantum SC fluctuation becomes essential, as in the example 
of Fig.l, due to a growth of A(0), one might expect the critical region of the thermal SC tran- 
sition in = to also widen due to the A(0)-growth. Actually, it was questioned in ref.26, 
through Figs. 7 and 8 there, whether the pictures arguing a pseudogap of SC fluctuation origin 
are consistent with the fact that the critical region of the H = transition does not widen 
much with underdoping. If the approach in §2 is extended to the H = case, however, we 
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find two theoretical facts satisfactorily explaining the sharp H = transition in underdoped 
cuprates. First, it is easily noticed that the microscopic T-dependence in the broad SC pseu- 
dogap regime of the coefficient of the gradient term in GL model results in a reduction of the 
width of the critical region around Tcq in H = 0. Actually, by noting that this GL coefficient 
in = is nothing but that of the i^-linear term in eq.(23), one easily finds that the 2D 
Ginzburg number nearTco is given, consistently with eq.(2), by 16Tr^[\{0)]'^ kB{Tco)^ /{4>odTQ), 
where the reduction factor (Tcq/Tq)^ arises from a microscopic T-dependence above Tco of 
the gradient term. Within the present approach in §2 where no specific origin leading to a 
reduction of the time scale 7 is assumed, this may be adequate for understanding the unex- 
pectedly narrow^^'^°^ H = critical region of underdoped cuprates. In a system with strong 
fluctuation, we have another mechanism of a shrinkage of the H = critical region due to the 
quantum SC fluctuation itself. An explanation of this mechanism will be given in Appendix 
B. This mechanism is particularly effective when the enhancement of quantum SC fluctuation 
is primarily due to a decrease of the dissipative time scale 7. In any case, a shrinkage of the 
H = critical region in a system with strong fluctuation is not a surprising phenomenon. 

We need to comment here on the definition of penetration depth comparable with exper- 
imental data.^^^ The actual penetration depth in = is defined as the mass of the gauge 
field through the gradient term. Then, by taking account of the T-dependence of the gradient 
term in the SC pseudogap regime mentioned above, the penetration depth to be observed near 
Tco is found to be not A(0)/(1 - T/T^if'^ but TcoA(0)/[To(l - T/Taif'\ If a T-dependence 
of the coefficient h in the SC pseudogap region can be neglected, the former is the quantity 
to be observed as the penetration depth in = and at low T where the fluctuation is 
negligible. Thus, in underdoped materials with wider SC pseudogap region, the [A(T)]~^ v.s. 
T curve is not linear even approximately, and, as observed in ref.59, the local slope \d\'''^ /dT\ 
increases on approacing Tco from below. However, we note that the A(0)-values we have es- 
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timatcd by fitting arc slightly longer than those estimated experimentally.^^^ If effects of the 
competing orders may be relatively negligible, the 6- value near Tco will be, as well as 70, larger 
than that near Tq. Thus, the differences between A(0)-values obtained through Fig. 6 and their 
experimental data^^^ do not necessarily require modification of the present theory. 

It should be stressed that key data-^^' showing the uncorrelated?^^ behavior between 

the Nernst coefficient (or the magnetization) and the resistivity were explained in the pre- 
ceding sections without taking account of electronic states peculiar to materials near a Mott 
transition. In the case of cuprates, such a behavior tends to arise in systems with low Hc{0) 
such as the overdoped and under doped materials, while it is rarely seen in optimally-doped 
cuprates with higher Hc{0)}^'^'^^ This nonmonotonic doping dependence in the cuprates we 
have explained in §3 strongly suggests that the uncorrelated behavior is not a consequence of 
a strong reduction of the friction coefficient due to a microscopic mechanism^*^'^^) peculiar to 
systems close to a Mott transition. It does not seem to us that the observation,^' -"^^'^^^ that 
the uncorrelated behavior and the sharp resistive drop are more remarkable in higher fields, 
while the fan-shaped broadening is usually seen in lower fields, can be explained in terms of 
such an approach^*^'^^^ based on a single vortex picture valid in low enough fields. Further, as 
emphasized in §1 and elsewhere, ^'^^ the uncorrelated resistive behavior is not peculiar to the 
cuprates near a Mott transition but also appears commonly in other superconductors such as 
the organic materials and disordered films with an ,s-wave pairing (see Figs. 9 and 13). 

In Fig. 6 (a), we have included a p{T) curve at 29(T) which is above H*2{0) (sec also 
Fig. 8). The sharp drop of this curve near 10(K) indicates a 3D VG transition point above 
H*2{T). Since, in principle, the VG transition may occur as far as the SC fluctuation is 
present, such a 3D VG transition and hence, a sharp resistive drop in the SC pseudogap region 
T*2{H) < T < Tq{H) may be present, in contrast to the scenario in ref.60 for the data,^) even 
in homogeneously disordered materials. The features seen in (K,Ba)Bi03 reported in ref.62 
may be a remarkable example with T*2 < T^g < Tq{H). 

It is interesting to connect the nonmonotonic^^^ doping dependence of T^{x) in LSCO to 
a sign change^^' of the fluctuation Hall conductivity a^^^xy It is now clear^^'^^^ that a sign 
reversal of Hall conductivity, usually seen in the vortex liquid regime, may occur above H*2{T), 
depending on the materials. It means that this Hall-sign reversal should be understood based 
on the fluctuation scenario^^^ unrelated to the electronic states of vortex cores. According to 
this scenario, the sign of a{[^xy is determined by that of dT^^^ /dx, and hence, if T^, is essentially 
identical with T^^, the sign change of dTjj/dx should appear directly as that of (7i[^xy, where 
was defined in §3. According to ref.63, the LSCO data in x > 0.12 show a Hall-sign 
reversal, while available p^y data in a; = 0.08^^^ have not shown any sign reversal, consistently 
with our expectation. On the other hand, according to ref.15, it appears that the Hall data 
of very underdoped Bi2201 still show a sign reversal, consistently with the monotonic doping 
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dependence of T^, of this material. This is why the identification between Ti, and T^^^ seems 
to be consistent with available Hall-resistance data. 

In hole-doped cuprates, the positive magnetoresistivity in the SC pseudogap region above 
Tco is enhanced with underdoping. Actually, this trend has led the authors in ref.47 to argue 
that ,^0 increases with underdoping. This enhanced magnetoresistivity above Tco can be seen 
as having a common origin to the resistivity curves^^' following the extrapolated normal 
curve even much below T*2{H) (< Tco). Actually, in a strong SC fluctuation case such that 
its quantum nature is no longer negligible, the Gaussian approximation for the SC fluctuation 
fails, and the interaction between the SC fluctuations may be important even much above Tcq. 
Hence, the quantum SC fluctuation may be a main origin of the enhanced^'') magnetoresistance 
in the SC pseudogap region (i.e., Tco < T < Tq). In any case, as demonstrated in rcf.20, 
focusing only on p(T) data in cases with low Hc{0) and hence, with large quantum fluctuation, 
tends to lead to an erroneous estimation of material parameters. A simultaneous study of other 
quantities, such as ?7^, measured consistently with p is indispensable. Further, the neglect of 
SC pseudogap region in thermodynamic quantities such as the magnetization^^-* has also 
erroneously led one to concluding a increasing with underdoping. 

During preparing this manuscript for submission, several related works ''^^ on transport 
phenomena in the vortex liquid regime were reported. The monotonic decrease^'^) of ^-c- 
companying underdoping we have concluded through fitting to data was also argued^'^^ from 
an extrapolation of Ncrnst data to very higher fields in which experimental measurements 
cannot be performed. The uncorrelated behavior between the resistivity and the Ncrnst data 
seen in underdopcd Pr2-xCcxCu04 was erroneously interpreted in rcf.74 as an inscnsitivity of 
the resistivity to SC fluctuations. As wc have clarified here, a theoretically valid explanation 
is provided only by noting the reduction of Gg induced by the quantum SC fluctuation at 
nonzero temperatures. 
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Appendix A 

Here, we give numerically useful expressions of the coefficients of the GL action in the 
case of s-wave dirty films.^^'^^'^^^ In obtaining them, we are largely based on the ordinary 
dirty limit in quasi 2D case where the 2D diffusion and fluctuation propagators are assumed 
consistently with a 3D electronic state (see the second paper of ref.23), and eflFects of Coulomb 
interaction between the electrons are included perturbatively and modelled in a form inter- 
polating between the low T/H and high T/H regions. The GL coefficients of the quadratic 
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term we use in §5 are given by 



2iy — 



(gi(0))-i~/xi-/xo (A-2) 
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where T^^ = O.UToH/H^^iO), Hf^{^) is Hq in dirty limit, i = T/T^^, R,, is the sheet 
resistance of a quasi 2D film at high T inversely proportional to the film thickness d, and 70 
in dirty limit with no electron repulsion was denoted as 7^'^). For simplicity, a SC pseudogap 
region is assumed to be absent (i.e., Tq = Tco), and 70 = 71 was additionally assumed because 
a difference between them is less important in dirty limit than that in clean limit (see eq.(23)) 
where 71 tends to vanish in T ^ limit. No interaction correction to ^i(O) was included 
because, as explained in §2 and ref.23, a detailed form of Gi{0) is not reflected in dfi and a^g. 

On the other hand, the coefficients b and bp, respectively, of the interaction and pinning 
terms will be expressed as 

' = Wr^5:(i + *> + V2))- (A-4) 

and 



V Rn / Stt 

^ ^ ^ n>Om>0 

X (1 + t{n + 1/2))-^ (1 + i{m + 1/2))-^ 

Effects of electron-repulsion on b and 70 were included altogether just in 70 since its effect 
on 70 was estimated in ref.53 to be much larger than that in 6, and they usually appears as 
the quantum fluctuation strength oc 6/70 in the quantum regime. We note that bp is Upfoo{0) 
in the notation of ref.23. Through the computations of resistivity curves shown in Figs. 12 
and 13, for simplicity, another parameter 27rTor, where r denotes the elastic scattering time, 
was fixed to 0.5. We expect a detailed value of this parameter not to significantly affect the 



32/?? 



J. Phys. Soc. Jpn. Full Paper 

numerical results. 

As a normal conductivity an for this case, the expression 

was used in the following figures, where Sq is given by cq.(lO). The second term of the nu- 
merator corresponds to the low T form of additional quantum SC fluctuation contributions^^^ 
to an excluded in the GL approach and was argued^^' '^^^ to be the origin of the fluctuation- 
induced'^^^ negative magnetoresistance in 2D and at low T. Such a behavior visible in higher 
H and lower T in Fig. 12 is a consequence of this fluctuation term. Further, to represent 
ln(ro/r)-term in the denominator of cr„ ensuring (T„(r 0) ^ 0, a form of an expected^^^ 
in the case with a strong spin-orbit scattering and with long-ranged Coulomb interaction was 
conveniently assumed as a model. In this way, when H and T are scaled, respectively, by 
-^c2(0) and To, both the GL coefficients and Rgdan are parametrized only by Rn/Rq- 



Appendix B 

Here, it will be explained how the quantum fluctuation itself reduces the width of critical 
region of the H = transition at Tc. To sketch the essence of this effect, we will just treat a 
counterpart of eq.(4) (with bp = AS; = 0) in the Hartree approximation and m H = 0. In the 
isotropic 3D case, the renormalized mass (corresponding to [^o(0)]~^) yields the relation 

/xfl - ln(r/rc) = J2[at^R;io) - C(0;a;)], (B-1) 

where 

4?>(r) = |16»2(A(0))V(f?&#Sl' is the 3D Gin.burg number a.t temperature T, and T, was 

redefined so that fiR{T = T^) = is satisfied. The thermal [lo = 0) part of the r.h.s. of eq.(Bl) 
becomes —{e"Q{T)Y/'^^iiR. If keeping only this contribution in the r.h.s. and focusing on the 
low ^iR limit, the critical exponent = 1 of the correlation length in the spherical limit is 
obtained. On the other hand, the sum of other a; 7^ terms in the r.h.s. is well approximated 

by —Tr~^fiR\J (7""^) if the quantum fiuctuation is strong enough. Then, solving eq.(Bl) 

(3) 

with respect to fiR, one finds the width of critical region to be estimated as Tce^ (Tc)/ [16(1 -|- 
TT~^\J Eq'^ (7^!))], implying that, with decreasing 7, the critical region is narrower. The basic 
reason of this narrowing of the critical region due to the quantum fluctuation consists in the 
fact that the width of quantum critical region at T = and in 3D is negligible, because the 
dimensionality of the dissipative quantum critical fluctuation at T = in L'-dimension is 
D + 2, which is above the upper critical dimension (i.e., four) for D > 27^^ 
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